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Abstract 



The existence of multi-pulse solutions near orbit-flip bifurcations of a primary single-humped 
pulse is shown in reversible, conservative, singularly perturbed vector fields. Similar to the non- 
singular case, the sign of a geometric condition that involves the first integral decides whether 
multi-pulses exist or not. The proof utilizes a combination of geometric singular perturbation theory 
and Lyapunov-Schmidt reduction through Lin's method. The motivation for considering orbit flips 
in singularly perturbed systems comes from the regularized short-pulse equation and the Ostrovsky 
equation, which both fit into this framework and are shown here to support multi-pulses. 
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1 Introduction 



The motivation for this work comes from the regularized short-pulse equation 

d z d t w + jw + d 2 w 3 + pdjw = 0, (1.1) 

derived in [TJ [12] as a model for short pulses in optical fibers, where w is the real component of the 
envelope of the electric field in the direction transverse to the direction of propagation, 7, (3 > are 
scaled real parameters, obtained from approximating the first order susceptibility tensor for wavelengths 
in the interval [1.6/um, 3.0/zm], and the Ostrovsky equation 

dzdtw + jw + d 2 z w 2 + pdjw = 0, (1.2) 

which was derived in (8] for nonlinear waves in the ocean in the small rotation limit and in [7j to 
describe oblique magneto-acoustic waves in rotating plasmas. More generally, the equation 

d z d t w + -fw + d 2 z w p + f3djw = 0, pSN, p > 2, (1.3) 



which encompasses (1.1) and (1.2), has been used in [6] as a model for nonlinear- wave phenomena in 



general rotating media. 



We are interested in travelling waves of the partial differential equation ( 1.3 ), which are solutions of the 



form w(z, t) = w(z — ct), where c is the wave speed. Using the travelling- wave coordinate y = z — ct, 



we find that travelling waves w(y) of (1.3) satisfy the ordinary differential equation 



{/3Wyy - CW + W P )yy + = 0. (1.4) 

It was shown in [T] that this equation supports single pulses (that is, homoclinic orbits) in appropriate 
parameter regions. Our goal is to show that this equation also admits travelling multi-pulses. 



The travelling-wave equation (1.4) associated with (1.3) is reversible under the operation y ^ —y and 



conservative. We will treat (3, which measures the strength of the regularizing fourth-order derivative 
term, as a small parameter. This turns the travelling-wave equation into a singularly perturbed problem 
and makes it amenable to a dynamical-systems analysis. In fact, we will see that the single pulses are 
in an or bit- flip configuration for ,9 = 0, that is, the single pulses are stronger localized than one would 
expect. For regular, not singularly perturbed problems, it is known that orbit flips often lead to multi- 
pulses [9] [TOl HU [13] . For regular, reversible, conservative systems, it was shown in [9] that a geometric 
condition decides whether multi-pulses exist: this condition basically measures whether the conserved 
quantity allows orbits to follow the primary homoclinic orbit near the origin within the energy surface 
of the homoclinic orbit. We will encounter a similar condition here for singularly perturbed equations. 



A key difference between (1.1) and (1.2), or more generally between even and odd values of p in (1.3), 
is equivariance under the transformation w \— ► —w, which holds only for odd values of p. In particular, 
the short-pulse equation has two single pulses, which are related by symmetry, while the Ostrovsky 
equation has only one single-pulse solution. Thus, we may expect that the short-pulse equation supports 
multi-pulses that follow the two single pulses in an arbitrary order as is the case for orbit-flips in regular 
reversible systems |11] . 
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To prove our existence and non-existence results, we combine results from geometric singularly per- 
turbation theory with Lyapunov-Schmidt reduction in the form of Lin's method. Lin's method is a 
tool for constructing multi-pulses that was introduced in [1] and further developed in |10J . In the 
travelling- wave context, it has been used, for instance, in [5] fT0 | HT] , IT5 ] . 

The rest of this paper is organized as follows. In fj2] we consider equation (1.4) in more depth and 
use it to motivate the slightly more general setup that we shall consider for our results. Our abstract 
multi-pulse existence results are stated and proved in ^3] and ^4] respectively. In §[5] we apply these 



results to (1.4). 



2 The profile equations 



As mentioned in the introduction, we treat the factor f3 in front of the fourth-order term in (1.4) as a 
singular perturbation parameter. This leads naturally to the scaling 



y, 



V = C i-PW, 



which turns (1.4), given by 



into the equation 



{fiWyy + W P - CW)yy + = 0, 



7 |/3| 

(sign^)^ + v p - sign(c)v) xx + —^v = 0. 

c 



Throughout this paper, we assume 

c>0, /3,7>0, 
so that the corresponding profile equation is given by 

(vxx + v p - v) xx + e 2 v = 0, 



1(3 



When e = 0, which corresponds to /3j = 0, equation (2.4) supports the pulse solution 



Z(*) 



P + 1\ P- 1 f p — 1 

sechp- 1 x 



We write (2.4) as the first-order system 



dx 



V-2 

V3 + V1- v p 

—evi 



in the coordinates 



(2.1) 
(2.2) 



(2.3) 



(2.4) 



(2.5) 



(2.6) 



V\ = V, V 2 = V , V3 = V + V P — V, 



1 When p is odd, our results apply also to the situation c,/3, 7 < 0, upon replacing v by —v. 
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The first-order system (2.6) is Hamiltonian with respect to the skew-symmetric operator 

/ -1 \ 
10 



J(e) 




e 







/ 



and the Hamiltonian 



H(v 1 ,v 2 ,v 3 ,v 4 ) := y 



V 



p+1 



— h viv 3 + — 

p+ 1 d 2 2 



Indeed, using V = (v i, v 2 , v 3 , v^) T , we have 



V = J(e)H v (V), 



where 



H V (V) 



( ^3 + Vi — ■! ' 

-v 2 

Vl 



? \ 



For e = 0, equation (2.6) admits the reversible homoclinic orbit 

Q(x) = (q(x),q'(x),0,0) T 



(2.7) 



with q(x) from (2.5 ). Note that this solution lies in the fast system of (2.6) and is therefore, by definition, 
in a singular orbit-flip configuration. It was shown in [1] by using geometric singular perturbation 
theory that this homoclinic orbit persists for e / 0. Furthermore, it was shown there that, for e ^ 0, 
the persisting homoclinic orbit lies no longer in the strong stable and strong unstable manifolds of 
the equilibrium V = 0, but converges to zero as |ac| — > oo with a smaller exponential rate that is 
proportional to e. 



3 Hypotheses and main result 



We consider singularly perturbed dynamical systems of the form 

= F(U,e), U = I : j G R z x R z , (3.1) 
where F(U, e) is a smooth nonlinearity, and (n, u) denote the slow and fast components of U. We now 




ef(u,u,e) 
g(u,u,e) 




list the assumptions on the system (3.1). 



Hypothesis 1 (i) The vector field is reversible, with reverser R(u\, 112, Ui, U2) = (^2, ui, U2, u\), so 
that F(RU, e) = —RF(U, e) for all (U, e). Furthermore, we assume that g(u, u, 0) = g(ui + u 2 , u). 

(ii) We assume that there exists a smooth function H(U,e) : M 4 x R — » R which is invariant under 
R and satisfies (Hu(U,e),F(U,e)) = for all (U,e) 6 I 4 x 1. We may normalize H so that 
H(0,e) = for all e. 



4 



Hypothesis [TJi) says that the system (3.1) is reversible: if U(x) is a solution of (3.1), so is RU(—x). 
Solutions with U(0) £ Fix(i?) satisfy U{x) = RU{—x) for all x and are referred to as reversible. 
Hypothesis [TJ^ii) means that (3.1) is conservative. A particular example of conservative systems are 
Hamiltonian systems U' = J(e)Hu(U,e), where J(e) : M 4 — ► M 4 is a skew-symmetric operator, such 
that J(e) = - J(e) T . We define 

£ e :={U; H(U,e) = 0} 



to be the zero level set of H. Our next assumption states that U = is an equilibrium of (3.1). 



Hypothesis 2 We assume that U = is an equilibrium of (3.1) for all e and that the linearization 
F;j(0, e) is given by 



Fu(0,e) 



where a(0) > 0. 



/ 


— e 


















e 


















-a{e) 







V 











a(e) 


) 



We remark that reversibility makes the spectrum of Fu(0,e) symmetric with respect to reflections 
across the imaginary axis. As a consequence of Hypothesis [2] and geometric singular perturbation 
theory [2], the system (3.1) has a two-dimensional center manifold near U = 0, and the flow on the 
center manifold is of the form 



u 




u G 



u + 0{\n\ 2 ) 

We focus on e > 0. In this parameter regime, there exists a unique smooth two-dimensional manifold 



VF u (0, e) of (3.1 ), which consists of the strong unstable foliation of the one-dimensional unstable man- 
ifold within the two-dimensional center manifold. For e > 0, this manifold coincides with the usual 
unstable manifold of U = 0. 

We will assume that the fast system, u' = g(0,u,0), has a homoclinic orbit Q(x,0) when e = 0. This 



orbit is automatically transversally constructed with respect to the full system (3.1 ), in the sense that 
W u (0, 0) and iy s (0,0) intersect transversally at (5(0,0) inside the level set Thus, the homoclinic 
orbit persists for e > 0. For e > 0, it may acquire a slow component u, and we shall assume this to be 
the case; see Figure [T] for an illustration. 



Hypothesis 3 We assume that (3.1) with e = has a reversible homoclinic orbit Q(x, 0) to the origin 



with Hjj(Q{0,0),0) 0. In this case, (3.1) has a reversible homoclinic orbit Q(x,e) to the origin for 
all e > close to zero, and we assume that 



lim Q(x,0)e x = e^, lim Q(x,e)e ex = rf{e)e\, 



?f(0) = 0, 



~de 



(0) + 0, 



where ej denote the canonical basis vectors. 
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Figure 1: The left panel illustrates the geometry of the center manifold and the unstable manifold 
W u (0, e). The right panel shows how the unfolding of the homoclinic orbit in the stable manifold for 
e > 0. 



In addition to the assumptions made above, the system (3.1 ) may be Z2-equivariant under the reflection 
U^-U. 

Hypothesis 4 We assume that F(U,e) is odd with respect to U so that F(—U,e) = —F(U,e) for all 
(U,e). 

If Hypothesis |4]is met, then Q and — Q are both homoclinic orbits, and we may seek iV-pulses that follow 
these two orbits in the order given by an arbitrary, but fixed, sequence {/%'}j=i,..jVj where Kj = ±1: 
the requirement is that the jth pulse in the TV-pulse follows kjQ; see Figure [2] for an illustration. We 
can now formulate our main result about the existence of multi-pulse solutions. 

Theorem 1 Suppose that Hypotheses^^and^are satisfied, and define 

a := sign {(Huu(0, 0)e 3 , Re 3 ) (Huu(0, 0)e 1 ,Re 1 }) . (3.2) 
If a = 1, then N -pulses with distances of order | lne| do not exist. If a = —1, then, for each N > 0, 



there exists an en > so that (3.1) has an N-pulse solution for each < e < e^v that winds N times 
around the primary pulse solution. The distances Lj between consecutive pulses in this N-pulse are 
given approximately by Lj w — lne as e — ► 0. 

Assume now that Hypothesis^ is also met. If a = 1, then N -pulses with distances of order \ In e| do not 
exist. If a = —1, then 2-pulses of up-up and up-down type exist for e > 0; see Figure^ Furthermore, 
for N > 2, each N-pulse either has k = ±(1, . . . , 1) or else at least one of the distances between 
consecutive pulses is not of order \ lne|. 

The condition a = — 1 was previously obtained in [9] for reversible, conservative orbit-flip bifurcations 



with hyperbolic equilibria. If the product of the scalar products in (3.2) is positive, then the energy 
inside the homoclinic loop in the fast and slow system has the same sign: this may prevent iV-pulses, 
which have zero energy, to pass through this area. 
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>0 JUUL 

/ (i) (ii) 




Figure 2: Panels (i) and (iv) contain phase portraits of homoclinic orbits in reversible systems with (iv) 
and without (i) 1%- symmetry. Panel (ii) illustrates the shape of a 3-pulse of up -up type with k = (1, 1, 1), 
while panel (Hi) contains an up-down-up 3-pulse with n = (1, —1, 1). 

4 Proof of Theorem [T] 



Throughout this section, we assume that (3.1) obeys Hypothesis El For a fixed sequence {^j}j 



--1,...N 

with Kj = ±1, we will then seek N- pulses for which the jth pulse in the iV-pulse follows KjQ. If 



Hypothesis [4] is not met, we simply set Kj = 1 for all j. 



4.1 Fenichel's normal form 



Hypothesis [2] implies that we can transform the singularly perturbed system (3.1) near the origin into 
Fenichel's normal form 



(4.1) 




/ e{A c {U c ,e)U c + B{U,e)[U ss ,U uu }) \ 




( UC \ 


A ss (U,e)U ss 


, u = 


U ss 1 £ R 2 x 1 x 


V A uu (U,e)U uu J 







where B(U,e) is a bilinear form for each (C7, e) and 

A c (0,e)= \-\ ° ), ^l ss (0,e) 




■1, 



A nn (0, e) 



1; 



see [2] and the references therein. Note that we rescaled x and e to normalize the fast eigenvalues to 
be ±1. We remark that the transformation can be chosen so that it is valid in the ball of radius 2 near 
the origin and respects reversibility and Z2 symmetry (when present). Since the slow center manifold, 
given by U ss = U uu 
manifold so that 

( A " iuc - e)u ' ), v=( ') e #, 

\ A u (U c ,e)U u I \ U u I 



- 0, is invariant under (|4.1|), we can straighten out its invariant stable and unstable 

a c (u c , e)u c = I - ) . / I i gpa 11.21 



where ^4 s (0,e) = — 1 and ^4 u (0,e) = 1 for all e. From now on, we will suppress the dependence of A J 
and B on the parameter e. 

We place two sections, Si n and S ou t, at f7 ss = 1 and U un = 1, respectively. We are interested in 
constructing solutions that need time Li for some large given Li to pass from S; n to S ou t- To find 
these solutions, we first construct convenient parameterizations of the two-dimensional stable and 



unstable manifolds of (4.1). We parametrize the two-dimensional stable manifold by Q + (x,b s ,e) so 



that Q + (L*, 6 s , e) S Si n for all (6 s , e) for some L* > 0, and Q + (x, 0, e) = Q(x, e) is the homoclinic orbit. 
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Furthermore, 6 s denotes the C/ s -component relative to the homoclinic orbit. Since the stable manifold 
has U u = U uu = 0, we easily find the expansion 

Q + (x, b s , e) = ([6 s + e + 0(|6 S | 2 + e 2 )]e~ ex , 0, [1 + 0(e) + 0{e- x )}e L *~ x , 0) , x > U 



from ( |4.l[ ) and (4.2): here, we also exploited Hypothesis |3| and assumed, without loss of generality, that 
dr/ s /de(0) = 1 (this can always be achieved upon replacing ej by — ej for j = 1,2 and rescaling the weak 
directions). The two-dimensional unstable manifold can be parametrized analogously by Q~(x,b,e), 
where Q~(— L*, 6 U , e) E Souti the scalar b u lies in the weak unstable direction U u , and the expansion 

Q~(x, 6 U , e) = (0, [6 U + e + 0(|6 U | 2 + e 2 )]e ex , 0, [1 + O(e) + 0(e x )]e L * +x ) , x < -L, 

holds. 

Fix a sequence {Ki}i=i,... N of numbers k« = ±1 that describes how a prospective 7V-homoclinic orbit 
follows the primary pulse Q(x,e) and its symmetric counterpart — Q(x, e). We shall also prescribe the 
times 2Lj that the individual pulses spend near the origin subject to the requirement that 

|e£;| < p 

for some sufficiently small <p< 1. We seek such iV-pulses using the parameterization 



(4.3) 



U l (x) = KiQ (x,bf,e) x e (— oo,0), 

ur( x ) = K i [Q-(x,bf,e) + V l -_ 1 (x)] ie(-Vi,0), 

U+(x) = Ki [Q+(x,ble) + V+(x)} xe(0,Li), 

U^(x) = K N Q + (x,b s N ,e) xG(0,oo) 

and require initially that 

U?(L.) e Kj S in , Ur(-L*) e K^ ont , U+(Li) = Ur +1 (-Li) 
for all i. These conditions mean that V^(x) should satisfy 

V+' SS (L*) = 0, Vp uu (-L*) = ( 4 - 4 ) 
K l [Q+{L u ble) + V l + {L l )} = Ki+1 [Q-(-L i ,bf +1 ,e) + Vr(-L i )}. (4.5) 

In addition, the functions K (x) should satisfy the system 

Vr x = F(Q-(x,ble) + V i Z 1 (x),e)-F(Q~(x,V,e),e), xG(-L,_ 1 ,0) 
V+ = F(Q + (x,ble) + V l + (x),e)-F(Q + (x,b!,e),e), xG(0,L 4 ), 

which we write as 

Vr x = F u (Q(x,b*,e),e)V- 1 + 0(\V- 1 \ 2 ) 1 xef-L^O) (4.6) 
V? = F u (Q(x,b s i ,e),e)V+ + 0(\V i + \ 2 ), x€(0,Lj. 

In deriving the above equations, we exploited the Z2-equivariance of the right-hand side F(U, e). 
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We focus on x > and solve the second equation in (4.6) for x £ (L#,Li). Using (4.1), we obtain the 
system 

K + ' C = e(A c (Q+> c (x, bl e))Vp c + 0{\V+> C \ 2 ) + B(V+)[Q + >%x, 6J, e) + V+' ss , ^ + ' uu ] 
V+' ss = A ss (Q + (x, bl e) + F+)(Q + ' ss (x, bl e) + V+' ss ) - A ss (Q + (x, b°, e))Q+' ss (x, 6f , e) (4.7) 
= A ss (Q + (x, bl, e))V+> ss + (A^(Q + (x, bl e))V+ + 0(| V+| 2 ))(U+' SS + Q + > ss (x, 6f , e))) 



We are interested in finding solutions V, (x) of (4.7) on (L*,Lj) so that 



(4.9) 



(v^ s ,v^ ss )(l*) = o, {vpW^U) = K,«D (4.8) 

for given small (a", a" u ). Proceeding as in f3j, we can construct these solutions and find that they obey 
the expansion 

V^' c (x) = e a i^afe 2 + 0(e< x - L ^\af\ 2 + €e-^ Li \af u \) 
V+> ss (x) = 0(e-^- s >\\V^ c \\+e-^ L '\ar\) 
V+> nu (x) = e a ™a? n , 

where \\V+' C \\ = sup xe[LmiLi] \vf' c (x)\ and 

af(x) = 6 f A c >»(Q + (y, 6|, e)) dy = e[x - U + 0(^(1^1 + e))] 
JLi 

(A™(Q+(y, bl e)) + 0(|^+(y)|)) dy = x-L l + 0(^(^1 + e)). 



«r(x) 

We also have 



11^11= OM + ee-^-^loTI 



(4.10) 



uniformly in Li. Proceeding in an analogous fashion for x < 0, we can solve the first equation in (4.7) 
with boundary data 



(yr>», VT' W )(-L*) = 0, (VP\ VP BB )(-Li) = (4, af) 



(4.11) 



and find that 



where 



V~' c (x) = e^a^ + Oie^+^lal 2 + ee-^- 5 ^\af 
Vr> uu (x) = 0(e (1 " <5)a: ||y-' c || +e- (1 - <5)Ll |af|) 
Vr' ss (x) = e a ^af, 

al(x) = e f X A c > s {Q + {y,V},e))&y = -e[x + L i + 0{L l {\V\+e))} 
J-Li 



(4.12) 



af(x) 
Again, we have 



(A ss (Q + (y, bl e)) + 0(|^+(y)|)) dy = -(x + L t ) + 0(L,(|a u | + e)). 



\\Vr> c \\ = 0(\al + ee-^ L >\ar\) 



(4.13) 



uniformly in Lj. 
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4.2 Construction of Lin orbits 

We now match the solutions Vi (x) first at x = L* in Si n and then at x = Li near the origin. To 
match in the section Ej n , we need to find expansions of the Poincare map that maps S ou t along the 
homoclinic orbit Q(x, e) into Si n . 

Lemma 4.1 The Poincare to Ei n satisfies 

n(U s , Q(-L*,e) + U u , U ss ; e) =: (0, Q(L*,e),0) + ft (4.14) 
= (0,Q(L„e),0) + [1 + 0(e)](U s ,U u ,U ss ) + 0((\U S \ + \U U \)\U\, (\U S \ + \U*\)\U, \U SS \ 2 ) 

and 

(Hu(Q(L*, e), e),ft) = eaU s + U ss + 0(e\U\(\U s \ + |C/ U |) + |C/ s f), (4.15) 
where we use the coordinates (U S ,U U ,U SS ) and (U s ,U u ,U un ) in S out and T,- m , respectively. 



Proof. The first expansion is really about the form of II when e = 0, since the e-dependent terms 



follow from Taylor expansion. For e = 0, equation (3.1) becomes 



u =0, u = g(u, u, 0). 
Hypotheses [T|i) an d [3] imply that the u-equation has a homoclinic orbit for all values of u, which 



implies that U ss = maps into U uu = for e = 0. The assertion (4.14) follows now from the slow-fast 



structure and Hypothesis [TJii) . Using the expansion of Q(x,e), we find that 

Hu{Q{L*,e),e) = eHuu(0, 0)ei + Huu(0, 0)e 3 = eae 2 + e 4 

Furthermore, since the gradient of the energy is perpendicular to the stable and unstable manifolds, 
which are parametrized by Q(—L*, b n , e) and Q(L*, b s , e), we see that there are no contributions of the 



form eb s or eb u to the scalar product (4.15). On account of the form of the error estimates for e = 0, 



the only nonlinear error terms that can appear in the expansion of the scalar product are as stated in 



(4.15). 



Evaluating (4.9) and (4.12) at x = dbL#, we obtain 



(1 + 0{eLi))afe 2 + 0{\af\ 2 + ee 



-(1-5)^, uu. 



-Li uu 



ani + 0(L 4 (K|+e))] 



and 



V t Z'i(-L*) 



(1 + 0(eLi_i))atiei + 0(|ati| 2 + ee' 



-(l-5)Lj_i I ss 



H-l\ 



-Li-i ss 



cdll + O^-idatil + e))], 



as well as V i ' (L # ) = V^_' 1 (— L*) = due to (4.8) and (4.11 ). To match the corresponding solutions, 



we transport Q (— L#, bf_ 1: e) + V r i _ 1 (— L*) from S out to Si n using the Poincare map n discussed in 
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Lemma 4.1 From (4.3), we find that the argument of II is given by 



Q-(-L*,bU,e) + V i Z 1 (-L*) 



( 



(l + 0(eL,_ 1 ))at 1 + 0(|at 1 | 2 + 6e- 



-(l-<5)Li_i I sb 



Q(-L*, e) + bf_ x + e + 0(|6ti| 2 + e 2 + |af_ x | 2 + ee" 

e -L ! _ la ^ i[1 + 0(Li _ i( | a s_ i | +e))] 



-(l-5)Li_i I ss 



H-l\ 



V 

We focus initially on the first two components of LT for which Lemma |4.1| gives 
n s ' u (Q-(-L*, bU, e) + Vr x {-L^ e) 



(4.16) 




(l + 0(eL,- 1 ))at 1 + 0(|af„ 1 | 2 + ee 
bf_ x + e + 0(|6til 2 + e 2 + K_il 2 + ee" 
+0 ((lat-il + l^-il + \e- L ^atMe + \aU\ + l^-il)) • 

anticipate the 8 
p\ hie | < Li < 



-(1— <5)£j_i I _ss 



**— II 



|„ss 



i-ll 



To simplify the calculations to follow, we anticipate the scalings we shall get: We shall choose 



which implies 



0(e), 



Q - Ll ss,uu 



0(6 2 ). 



(4.17) 
(4.18) 



Using these estimates, which we shall verify later in the proof, equation (|4.16|) becomes 
lT> u (Q-(-L*,^i,e) + Vr 1 (-L*);e) = 



+ 



(l + o(l))af_ 1 + 0(e 2 ) 
^-i + 0(e 2 ) 





^ Q(L*,e) 

Expanding the first two components of Q + (L*, bf, e) + VJ + (L*), and setting them equal to the compo- 
nents of (4.16), we arrive at the equations 

b\ + 0{e 2 ) = (l + o(l)K_ 1 + 0(e 2 ) 
&"-i + 0(e 2 ) = (l+o(l)K + 0(e 2 ), 

which we can solve by the implicit function theorem to get 

b\ = (1 + o(l)K_i + 0(e 2 ), &VLi = (1 + o(l))< + 0(e 2 ). 

It remains to match the third components in the ?7 ss -direction in Sj n . Projecting the difference of 
£+(£*,&?, e) + F+(L*) and U S > U (Q~~(-L*, 6£_ l3 e) + V^Z^-L*); e) onto Hu(Q(L*,e),e), we obtain 

& := (H u (Q(L„e),e),U s ' u (Q-(-L„bt 1 ,e) + Vr^-LJ; e) - Q+(L«, 6?, e) + F+(L*)> (4-19) 



e ( r(l + o(l))K-at 1 ] + (l + o(l))[e 



-Li_i,,ss 



afii]+0(e J 



We will solve the equations = at the very end of our analysis. 



Next, we match the piecewise defined solutions at x = Li. Evaluating (4.9) at x = Li and using (4.10), 
we obtain 



V^iLi) 
V^(L t ) 



a^ + 0(e 2 ) 
0(e 2 ^) 
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Similarly (4.12) and (4.13) imply 



Vr' C (-Li) = a^ + 0(e 2 ) 

V-> ss (-U) = af 

v-> uu (-U) = o( e 2 - 5 ). 

The matching condition ( |4.5| ), 

Ki[Q + {L h b% e) + VfiU)] = K i+1 [Q~(-L u bf +1 ,e) + Vr{-Li)], 

then becomes 

/ (e + ^ + 0(e 2 ))e- eL " \ / o| + 0(e 2 ) \ 



ft.. 



< + 0(6 2 ) 

e" L * +0(e 2 - 5 ) 



We find 



and 



(e + 6- +1 + 0(e 2 ))e-^ 
af 

V «r y v e-^ + o(e 2 - 5 ) y 

af = ft^+ie"^ + 0(e 2 ^), a? u = ft^ m e- L < + 0(e 2 " 5 ) 



a? = KiK i+1 (e + ati + 0(e 2 ))(l + o(l)), af = Ki K i+ i(e + af +1 + 0(e 2 ))(l + o(l)). 

We can solve the preceding system of equations to get 

/ i \ I N \ 

a\ = en i+ i o(l) +y\j , af = eKi o(l) + K j ' i = l,...N — 1. 

\ j=l J \ j=i+l I 



Note that this validates the first relation in (4.18) that we used above. We now substitute these 



expressions into the remaining bifurcation equations (4.19), which become 



6 = 6a(l+o(l))K-at 1 ] + (l + o(l))[e~^ar-e-^- 1 a- 1 ]+0( e J 

N i-1 



(4.20) 



e 2 K i( r(l+o(l)) 



for i = 1, . . . , JV, where Lq = Ljv = 00. 



+ (1 + o(1))k; [K m e" 2L * - Ki-ie" 2 ^- 1 ] + 0(e 3 ) 



4.3 Bifurcation equations 



On account of Lemma 3.2], it suffices to solve (4.20), given by 

i-1 



eV(l + o(l)) 



N 



j=i+l j=l 



+ (1 + o(l)) [K l+1 e~ 2U - Ki-ie- 2 ^- 1 ] + 0(e 3 ) = 0, 



where i = 1, . . . , N — 1, since £j\r then vanishes automatically due to the presence of the conserved 
quantity H that we assumed to exist. Define > via 



e 2 a ?; = e- 2Li , 



12 



and observe that this yields (4.17) and the remaining second estimate in (4.18) that we used to derive 



the bifurcation equations. We obtain 

N 



*(1 + 0(1)) 



i-1 



j=i+i j=i 



+ (1 + o(l)) [K i+1 ai - Ki_iOi_i] + 0(e) = 0, 



(4.21) 



where i = 1, . . . , N — 1. Setting e = 0, we arrive at the system 



Ki+idi — Ki-iai-i — a 



i-1 



N 



j=l j=i+l 



(4.22) 



or, written out in detail, at 



K 2 a\ 

— Ki-lCLi-1 + Ki+ldi 
-K7V_ 2 aAT_2 + KNdN-1 



-<t(k 2 + . . . + Kjv) 

-a(n i+ i + . . . + k n ) + cr(«j_i + . . . + «i) 

-(JKN + a{KN-2 + . . . + Kl) 



for i = 2, . . . , N — 2. The left-hand side is linear and invertible in a = (ai, . . . , ajv_i), as it corresponds 
to a lower triangular matrix with strictly positive entries Ki+i on the diagonal. Thus, if we can find 
a positive solution a* > to (4.22) for a given sequence Ki, then we can solve the full equation (4.21) 
using the implicit function theorem. 



First, we shall look for TV-pulses of up-up type and therefore set Ki = 1 for all i, so that (4.22) becomes 

at - Oi-i = -a(N + 1 - 2i), i = 1, . . . , N - 1. 

The equation for i = 1 is a± = —a(N — 1), which has a positive solution only when a = —1. Thus, 
iV- pulses of up-up type can exist only for a = — 1. Hence, we take a = — 1 and therefore need to solve 

Oi - di-i = N + 1 - 2i, i = 1, ...,N- 1. 

This system has the positive solution = i(iV — i) > with i = 1, . . . , N — 1, since 

m - cii-i = i(JV -»)-(*- - (i - 1)) = i(N - i) - (i - l)(iV - i + 1) = JV + 1 - 2». 

The preceding discussion therefore shows that there is a iV-pulse of up-up type for each e > and that 
there are no other iV-pulses of up-up type whose distances satisfy (4.17). 

Next, we investigate 2- pulses, when Hypothesis [4] is met. In this case we need to solve the single 
equation 

K2d\ = — aK 2 

which has the solution a± = —a, independently of k 2 - Thus, for a = — 1, 2-pulses of both up-up 
{k\ = k 2 ) and up-down {k\ = —k 2 ) type exist. 

It remains to consider the case N > 2 for arbitrary sequences {/tj}. 
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Lemma 4.2 Without loss of generality, let K2 = 1, and assume that m E {±1} for i = 1, . . . , N, where 
N > 2. // (4-22) has a solution a = {ai\i=\ : ,,^N-i with a» > for all i, then necessarily a = — 1 and 
Ki = 1 for all i. 



Note that this completes the proof of Theorem [T] due to the arguments presented above, upon taking 



the assumption (4.17) about the distances between consecutive pulses into account. 



Proof. Let 



i-1 



N 



bi = a 

Using induction, we find that 

h = h + a 



j = l j=i+l 



i = 1,...,N- 1. 



(4.23) 



ki + Kj + 2 , 



i=2 



i > 2. 



Similarly, an induction argument shows that the solution of (4.22) is given by 



i = l,...,JV-l. 



(4.24) 



Using the normalization k<i = 1, we find that a\ = K2^i = fti > 0, hence &i > 1. We now proceed again 
by induction and make the following induction statement at the ith step: 



J^Kjbj =i[bi + a(i-l)}, 



h > 



1 when (j = l 
i when a = — 1, 



Kj = 1 for j < i + 1. 



Using (4.24) and the assumption that aj > for all j, it is not difficult to check that the statement is 



true for i = 2 and to carry out the induction step from i to i + 1, and we therefore omit the details. 
Thus, we find that Kj = 1 for all i, which implies b\ = —a{N — 1). When used in combination with 
b\ > 1, these statements prove the lemma. ■ 



5 Application to the generalized Ostrovsky equation 



We now return to the generalized Ostrovsky equation (1.3) 



d z dtw + + d 2 z q p + Pdjq = 0, p G N, p > 2. (5.1) 
As we had seen in ^2j travelling waves of the form w(z, t) = w(z — ct) satisfy the differential equation 



(2.2) 



(fiVyy +V P ~ CV)yy +7^ = 0, 



where y = z — ct is the travelling coordinate. The scaling (2.1) and the condition (2.3), namely c > 



and P, 7 > 0, allowed us to write this equation as (2.4) 



{vxx + v p - v) xx + e v = 0, 



2 n 2 7/3 
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which supports the pulse 



p + 1 \ p- 1 ( p — 1 
sechp- 1 x 



(5.2) 



when e = 0. Equivalently, we can write this equation as the first-order system (2.6) 





( 




\ 


/ 


V2 




d 




v 2 






V3 + V1- 


< 


dx 










—ev^ 






V 


V 1 


/ 


V 


ev\ 





(5.3) 



Equation (5.3) is reversible with respect to the reverser R(v±, V2, i>3, W4) = {v\, — V2,vs,— V4) and con- 



servative with respect to the first integral (2.7), 

«? 

H(Vl,V2,V 3 ,V4) 



V 



p+1 



V U1U3 1 . 

2 p+1 2 2 



and it is now easy to check that Hypothesis [T] is met. It was shown in [T] that Hypotheses [2] and [3] are 



met for the pulse given in (5.2). Furthermore, we have 



D 2 H(0) 



( 


1 





1 





\ 







-1 












1 













V 











1 


/ 



and, using that the slow and fast stable eigenvectors of the linearization of (5.3) about the origin are 
given by 



(1,0,-1, -if, 



(1,-1,0, of , 



we find that 



sign((D 2 iJ(0)e s ,i?e s ) (B 2 H{0)e ss , Re ss )) = -1. 

—V if and only if p is odd. We can therefore apply 



Finally, ( |5.3[ ) is equivariant with respect to V 
Theorem [T] to get the following result: 



Theorem 2 Fix 7, c > 0. For each N > 2, there exists a (5m > such that the generalized Ostrovsky 



equation (5.1) has an N -pulse wn(z — ct) of up-up type for each (3 with < /3 < /3jv- The distances 
between consecutive pulses in the N -pulse wn(z — ct) are of the order — /31n/3. If p is odd, —wn{z — ct) 
is also an N -pulse, and 2-pulses of both up-up and up-down type exist for < j3 < fo- 
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